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Negatively charged hadron (h~) production in Au+Au collisions at BNL-RHIC is stud- 
ied with the STAR experiment. Results are presented on h~ multiplicity, pseudorapidity 
and transverse momentum distributions at ^/snn = 130 GeV. 

1. Introduction 

Global observables such as the multiplicity and the inclusive single particle transverse 
momentum (p±) and pseudorapidity (77) distributions of hadrons have been valuable tools 
in studying heavy ion collisions. They represent the system at kinetic freeze-out, late in 
the evolution of the system when hadrons no longer interact with each other, and their 
momentum spectra do not change further. These final-state observables supply essential 
constraints on the possible evolutionary paths of the system that can help establish condi- 
tions in the early, hot and dense phase of the collision. A discussion on global observables 
can be found in . We summarize here results on the minimum-bias multiplicity, p± and 
rj spectra of negatively charged hadrons. 

2. Experiment and Analysis 

The STAR experimental setup for the first RHIC run is described in ||. The analysis 
is based on charged particle tracking in the STAR Time Projection Chamber (TPC). 
Triggering was achieved using two hadronic calorimeters (ZDCs) in the very forward 
region and an array of scintillator slats arranged in a barrel (CTB) around the TPC. 

Particle production was studied through the yield of primary negative hadrons (h~), 
i.e. 7r~, K~ and p including the products of strong and electromagnetic decays. Low 
momentum (< lGeV/c) particle identification was done via dE/dx in the gas of the TPC. 
Only negatively charged hadrons were studied in order to exclude effects due to participant 
nucleons. Charged particle tracks reconstructed in the TPC were accepted if they fulfilled 
requirements on number of points on the track and on pointing accuracy to the event 
vertex. The measured raw distributions were corrected for acceptance, reconstruction 
efficiency, contamination due to interactions in material, misidentified non-hadrons, the 
products of weak decays, and track splitting and merging. The tracking efficiency was 
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found by embedding simulated tracks into real events at the raw data level, reconstructing 
the full events, and comparing the simulated input to the reconstructed output. 

For this analysis, the acceptance for tracks within the fiducial volume having p± greater 
than 0.3 GeV/c is found to be 95%. The tracking efficiency is found to be between 70-95%, 
depending on the particle p± and the total multiplicity of the event. 

3. Results and Discussion 

Figure [l] shows the h~ multiplicity distribu- 
tion for minimum-bias Au+Au collisions with 
\r]\ < 0.5 and p± > 100 MeV/c. The data 
were normalized assuming a total hadronic 
inelastic cross section of 7.2 barn for Au+Au 
collisions at y/s NN = 130 GeV, derived from 
Glauber model calculations ||[| . 

The systematic error on the vertical scale is 
estimated to be 10% and is dominated by un- 
certainties of the total hadronic cross-section 
and the shape at low N h - . The systematic er- 
ror of 6% on the horizontal scale is depicted 
by horizontal error bars on a few data points. 
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by the collision geometry: large cross section 
Figure 1: Multiplicity distribution of h~, at low multiplicity, corresponding to large im- 
shaded area is for the 5% most central col- pac t parameter, followed by a region of slowly 
lisions. falling cross section over a wide range of mul- 

tiplicity as the nuclei overlap, and a rapid decrease for near head-on collisions, where 
the shape is determined by fluctuations in the collision geometry and particle production 
process given the finite detector acceptance. The distribution for the 5% most central col- 
lisions (360 mbarn), defined via ZDC coincidence, is shown as the shaded area in Fig. [TJ. 

Figure 0, upper panel, shows the h~ p± distribution for the 5% most central colli- 
sions at midrapidity (|?7| < 0.1). Only statistical errors are shown; the systematic er- 
rors are estimated to be below 6%. The data are fit by a power-law function of the 
form (1/px) dN h -/dp± = A (1 + p±/po)~ n - As reference data sets, we also show the p±- 
distribution of negatively charged hadrons for central Pb+Pb collisions at y/s NN = 17 GeV 
from NA49 || and for minimum-bias pp collisions at yfs= 200 GeV from UAl |J, fitted 
with the same function. The UAl invariant cross-section Ed 3 a/d 3 p reported in Ref. || 
was scaled by 27r/<7 ine i, where <r ine i = 42 mb. The p± spectrum measured by STAR reflects 
a systematic increase in (p±) compared to that from both central A+A collisions at much 
lower energy (NA49) and pp collisions at a comparable energy (UAl). 

The lower panel of Fig. |2| shows the ratio of the STAR and UAl ^-distributions. Since 
the UAl distribution is measured at y/s = 200 GeV, the invariant cross-section in each 
p_i_-bin is scaled by two factors for a quantitative comparison to the STAR data: (i) 
i?(130/200), the p^-dependent ratio of the h~ yields in pp collisions at \fs = 130 and 200 
GeV, and (ii) Taa= 26 mb" 1 , the nuclear overlap integral |§ for the 5% most central 
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Au+Au collisions. See |7j for further details. 

There are two simple predictions for the 
scaled ratio. The study of lower energy col- 
lisions has shown that the total pion yield 
due to soft (low p±) processes scales as the 
number of participants, i.e. , "wounded" 
nucleons (iV WN ) in the collision, (e.g. 0). 
The scaled ratio in this case is 0.164 as- 
suming 172 participant pairs |4]] and a mean 
number of binary collisions (iVec) of iVec = 
o"inei Taa = 1050 for the 5% most central 
Au+Au events. We assume <7 ine i = 40.5 
mb for pp collisions at yfs = 130 GeV. If 
hadron production is due to hard (high p±) 
processes and there are no nuclear-specific 
effects (see below), the hadron yield will 
scale as Nbc, which is proportional to the 
nuclear overlap integral Taa- In this case 
the value of the ratio is unity. There are 
important nuclear effects which should al- 
ter the scaling from these simple predic- 
tions, including e.g. initial state multi- 
ple scattering |J, jet quenching |l(J, and 
radial flow 
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Figure 2: Upper panel: h~ p_i_-spectra for the 
5% most central collisions. NA49 central and 
UA1 data are also shown. Lower panel: ratio 
of STAR and scaled UAl ^-distributions. 



Each of these exhibits 
characteristic features as a function of p± 
and system size. The scaled ratio shows 
a strong p± dependence, starting close to 
iVwN scaling at low p± and approaching 
iVec scaling, not reaching the latter in the 
measured p± range even including the er- 
rors from the _R(130/200) and Taa scaling 
shown in gray. This behaviour is consistent with the presence of radial flow, as well as 
the onset of hard scattering contributions and initial state multiple scattering with rising 
p±. The shape of the distribution at higher p± is discussed in 0. 

The h~ density at midrapidity for p± > 100 MeV/c and \r]\ < 0.1 is dN/di]\ v=Q = 253 ± 
l s tat. ± 16 S2/st .. Extrapolation to p ± = yields dN/ dr]\ n=0 = 275 ± l stat . ± I8 syst ,. Assuming 
an average of 172 participant pairs per central Au+Au collision, this corresponds to 
1.60 ± 0.13 h~ per participant nucleon pair per unit of pseudorapidity. This is a 35% 
increase over pp collisions at the same energy [l^| . Comparison with central Pb + Pb 
collisions at SPS shows an increase in the h~ yield per participant of ~ 49% at RHIC. 

Fig. H shows the centrality dependence of the r\ distribution. One sees the expected rise 
in particle yield with collision centrality. Preliminary analysis shows very little difference 
in the shape of the distribution at midrapidity with increasing centrality. Using again a 
power-law fit to the p± spectra we obtain (p±) which is shown in Fig. |asa function of 
centrality. Errors on the horizontal scale are systematic. We see an increase in (p±) of 
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15% from the most peripheral to the most central events. The values for the UAl and 
NA49 reference data sets are also shown. 
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Figure 3. h~ 
centralities. 



rj distribution for different Figure 4. Centrality dependence of (p±). 



4. Conclusion 

We have presented preliminary results on charged particle production at RHIC. We 
find dNh-fdrj per participant in central Au + Au collisions at ^/s NN =130 GeV increases 
by 35% relative to pp and 49% compared to nuclear collisions at y/s NN = 17 GeV. The p± 
distribution is harder than that of the reference systems for the p± region up to 2 GeV/c, 
with the most central events having the highest (p±). Scaling of produced particle yield 
with number of participants shows a strong dependence on p±, with Wounded Nucleon 
scaling achieved only at the lowest measured p±. 
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